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I.  INTRODUCTION 


To  properly  describe  chemical  reactions 
in  non-permeable  porous  media,  the  "pore  tree” 
was  introduced  by  Simons  and  Finson  (1979) 
and  Simons  (1982).  This  pore  structure  was 
developed  via  analogy  to  the  kinetic  theory  of 
gases:  the  pore  length  is  analogous  to  the  mean 
free  path  and  under  the  assumption  that  the 
pore  aspect  ratio  (pore  length  to  radius)  is  a 
constant,  a  pore  size  distribution  was  obtained 
that  has  been  confirmed  for  coal,  coal  char, 
sorbents,  catalysts  and  even  for  kidney  stones 
from  both  men  (low  porosity  oxalate)  and 
women  (high  porosity  phosphate).  The  pore 
tree  was  statistically  derived  from  the  pore  size 
distribution  and  allows  the  orderly  migration  of 
a  reactant  gas  from  the  large  pores  to  the  small 
pores  (see  Fig.  1).  A  detailed  description  of 
the  pore  tree  and  the  coupled  transport  and 
chemistry  is  given  by  Simons  (1982,  1983a). 

The  spatially  dependent  transport/reaction  equations  are  solved  for  a  single  pore  tree  and  then 
the  total  contribution  of  all  trees  (of  all  sizes)  in  the  system  is  obtained  by  summing  the 
contribution  of  each  tree  that  reaches  the  exterior  of  the  system.  This  is  distinct  from  the 
"bulk"  transport  equations  that  integrate  over  all  pores  at  a  fixed  point  in  space  before  they  . 
integrate  spatially.  The  "bulk"  transport  approach  inherently  assumes  instantaneous  mixing 
between  large  and  small  pores  at  a  single  point  in  space  (Random  Pore  Models)  while  the 
pore  tree  relaxes  this  assumption  with  a  structure  that  is  developed  statistically  and  is 
relatively  easy  to  work  with  mathematically. 


1 

mtm  ■ 

c 

wwwwwwww^ 

rp  =  Appropriate 
Mean  Radius  at  x 


Ummmm  Particle  Surface 


Figure  1.  The  Pore  Tree 


In  all  previous  applications  (sorbents,  catalysts  and  coal  combustion)  the  pore  tree 
represents  an  isolated  sub-structure,  allowing  diffusion  into  and  out  of  the  porous  media 
without  permitting  convection  through  the  media.  Under  the  current  ARO  program,  Simons 
Research  Associates  is  adapting  its  pore  tree  model  to  describe  the  permeable  pore  structure 
and  convective  transport  characteristic  of  the  subsurface  flow  of  water  in  soil,  the  dispersion 
of  contaminants  and  the  in-situ  remediation  of  contaminated  sites.  The  random  nature  of  the 
pore  structure,  which  formed  the  basis  of  the  statistical  derivation  of  the  pore  tree,  is  applied 
to  the  underground  structure  of  porous  soil,  sand  and  rocks.  The  statistical  description  of  the 
pore  tree  is  extended  to  describe  the  concept  of  permeability,  thus  allowing  convective  flow 
as  a  transport  mechanism  in  addition  to  the  diffusive  transport  already  considered.  Physically, 
this  amounts  to  statistically  determining  the  "branches"  that  are  common  to  several  trees  to 
allow  percolation  through  the  large  scale  (mobile)  structure  in  addition  to  diffusion  through 
the  smaller  scale  (immobile)  structure.  While  permeability  is  dominated  by  the  largest  pores, 
it  is  important  to  determine  the  level  of  convection  that  is  occurring  at  the  intermediate  scales 
in  order  to  accurately  describe  transport.  The  resulting  pore  structure  will  provide  an  analytic 
description  of  an  underground  network  upon  which  transport  and  coupled  chemical  reactions 
may  be  accurately  superimposed. 
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II.  ISOLATED  PORE  TREE:  THE  STRUCTURE 

Following  the  pore  structure  theory  of  Simons  and  Finson  (1979)  and  Simons  (1982), 
consider  a  spherical  porous  particle  of  radius  a,  containing  pores  of  length  Ip  and  radius  rp. 
The  pore  dimensions  range  from  a  microscale  of  the  order  of  Angstroms  to  a  macroscale 
which  is  a  significant  fraction  of  the  particle  radius.  The  radius  of  the  largest  pore  is  denoted 
by  r^  and  is  given  by 


where  0  is  the  total  porosity  of  the  particle  and  K0  is  a  constant  of  integration, 
approximately  equal  to  five,  which  relates  the  pore  length  to  its  radius 

In  (2) 

The  radius  of  the  smallest  pore  is  denoted  by  r^  and  is  given  by 


^“20/pp^ 

where  p3  is  the  density  of  the  solid  matrix,  Sp  is  the  specific  internal  surface  area  (several 
hundred  m2/g),  and 

P=ln(rJrJ 


The  particle  contains  a  continuous  distribution  of  pore  sizes  from  rmin  to  rimx  .  The 
number  of  pores  within  an  arbitrary  plane  of  cross-sectional  area  A  and  with  radius  between 
rp  and  rp  +drp  is  denoted  by  g{r^Adrp .  The  pore  distribution  function  g(/p  is  given  by 


5(rp=e/2*pr,3  (5) 

where  g(rp)  indicates  an  average  over  all  inclination  angles  between  the  axis  of  the  pore  and 
the  normal  to  the  plane.  Due  to  the  random  orientation  of  the  pores,  the  intersection  of  a 
circular  cylinder  with  a  plane  is  an  ellipse  of  average  area  2jrrp2.  Hence,  the  porosity  is  the 
2rcrp2  moment  of  g{r^  and  the  internal  surface  area  is  the  47trp  moment  of  g{r^ .  The 
expression  for  g{r  )  was  derived  (Simons  and  Finson,  1979)  from  statistical  arguments  and  has 
been  validated  through  extensive  comparison  of  the  predicted  volume  and  surface  area 
distributions  with  mercury  intrusion  data  (Stacy  and  Walker,  1972).  This  has  been 
accomplished  for  coal,  char  derived  from  that  coal  (Kothandaraman  et.al.,  1984),  sorbents, 
catalysts  and  even  kidney  stones  from  both  men  (oxalate)  and  women  (phosphate). 
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A  characteristic  feature  of  the  l/rp3  distribution  depicts  that  the  pore  volume  between 
r^  and  rp  increases  linearly  with  the  natural  log  of  rp.  It  is  the  functional  form  of  this 
relationship, 


Pore  Volume 


r 


In/ 


(6) 


that  depicts  the  inverse  cubic  dependence  of  g(r p  on  rp.  A  linear  display  of  mercury  intrusion 
volume  vs.  ln(rp)  always  infers  a  l/rp3  distribution. 

The  number  of  pores  within  the  bulk  volume  V  whose  pore  radius  is  between  rp  and 
rp+drp  may  be  defined  by  Vfir^drp.  The  pore  volume  is  expressed  as  the  Ttrp2lp  moment  of 
firj  and  the  internal  surface  area  is  the  27crplp  moment  of  fir^ .  The  pore  size  distribution 
functions  (  fir^  and  g(rp  )  are  clearly  not  independent.  The  definitions  of  porosity  and 
internal  surface  area  infer  that  fir p  is  related  to  g(r^  by 


gir^ArJip 


(7) 


Equation  (7)  simply  states  that  the  probable  number  of  pores  intersecting  an  arbitrary  plane 
increases  with  the  length  of  the  pore  and  with  the  density  of  pores. 


The  length  of  a  pore  is  determined  by  an  arbitrary  intersection  with  another  pore  and 
is  expressed  (Simons  and  Finson,  1979)  as  a  collision  integral  over  the  pore  distribution 
functions.  The  analysis  suggests  that  lp  ,  g(r^  and  fir p  are  proportional  to  rp  ,  l/rp3  and  l/rp4 
respectively.  The  constants  of  proportionality  are  obtained  from  integral  constraints,  i.e.,  the 
total  porosity  and  internal  surface  area  contained  in  the  pore  structure.  The  expression  for 

firj  is  given  by 


(8) 


where  the  constants  were  defined  above. 

The  pore  volume  distribution  corresponding  to  these  distribution  functions  is  similar  to 
that  utilized  in  the  random  pore  model  (Gavalas,  1980  &  1981).  However,  the  pore  tree 
model  and  the  random  pore  model  differ  dramatically  in  their  choice  of  the  pore  aspect  ratio 
(length  to  diameter)  and  its  implications  with  respect  to  pore  branching.  The  random  pore 
model  allows  a  single  pore  to  connect  two  larger  pores.  This  picture  lends  itself  to  the 
idealization  of  instantaneous  mixing  between  the  pores  and  requires  that  the  pore  aspect  ratio 
be  of  the  order  of  one  hundred.  The  pore  tree  theory  uses  data  for  r^  to  imply  (via  K0)  that 
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all  pores  possess  an  aspect  ratio  of  the  order  of  ten.  Hence,  small  pores  may  connect  to  larger 
pores  only  on  one  end  and  all  pores  must  branch  from  successively  larger  pores  like  a  tree  or 
river  system. 

Each  pore  that  reaches  the  exterior  surface  of  the  particle  is  depicted  as  the  trunk  of  a 
tree.  The  size  distribution  of  tree  trunks  on  the  exterior  surface  of  the  particle  is  denoted  by 

g{r^4na2drt  where  g(rt)  is  functionally  identical  to  g(r^ .  Each  trunk  of  radius  rt  is 
associated  with  a  specific  tree-like  structure.  Let  Nt  be  defined  as  the  branch  distribution 
function  where  Ntdrp  is  the  number  of  pores  of  radius  rp  (within  size  range  drp)  in  a  tree 
whose  trunk  radius  is  rt .  The  total  number  of  pores  of  radius  rp  in  a  sphere  of  radius  a 

may  be  expressed  as  4/3 na2f(r^drp  or,  as  the  sum  of  all  pores  of  radius  rp  contained  within 
every  tree  in  the  porous  sample,  plus  all  pores  of  radius  rp  that  are  themselves  the  trunk  of  a 
tree.  Hence, 


-^7t£j3/rp)=  f  Ng(rt)4na2drt  +  4Tia2g(rp) 


(9) 


where  g(r)  is  the  number  of  tree  trunks  per  unit  external  area  of  the  porous  sample  and  only 
those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp .  Using  the 
previously  derived  expressions  for  r^,  g(fp)  and  ,  Eq.(9)  is  identically  satisfied  by 


xi  3/  * 
N,=rt/rp 


(10) 


The  branch  distribution  function  completely  characterizes  the  pore  tree.  The  internal 
surface  area  and  pore  volume  associated  with  each  pore  tree  are  denoted  by  St(r,)  and  Vt(r,), 
respectively,  and  are  expressed  as  the  sum  of  the  contributions  from  the  trunk  and  that  from 
the  branches. 


St(rt)=2nr/(+ f  2xrplpN/trp 


(11) 


*mlii 


V'(rt)=nr?l  + f  nr2pl^/irp 


(12) 


Using  Eq.(10)  for  Nt,  S/rJ  and  Vt(r,)  become 
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(13) 


W=2  *rh 


(14) 


where  the  (1-0)  term  in  St  has  been  included  to  account  for  pore  combination  (Simons, 
1979a). 


The  surface  area  associated  with  the  pore  tree  may  be  several  orders  of  magnitude 
greater  than  the  surface  area  of  the  trunk.  However,  the  volume  of  the  pore  tree  may,  at 
most,  be  one  order  of  magnitude  greater  than  that  of  the  trunk.  It  should  also  be  noted  that 
the  above  expressions  for  St  and  Vt  reduce  to  those  appropriate  to  a  single  cylindrical  pore  in 
the  limit  of  rt  ->rmin  (the  leaf  of  the  tree).  Furthermore,  the  integrals  of  St(r,)  and  Vt(r,)  over 

all  g(r)  recover  the  total  internal  surface  area  and  pore  volume  of  the  porous  sample. 

Each  trunk  of  radius  rt  is  associated  with  a  specific  tree-like  structure  with  continuous 
branching  to  ever  decreasing  pore  radii.  The  radius  and  number  of  pores  is  a  unique  function 
of  the  distance  x  into  the  tree.  The  coordinate  x  is  skewed  in  that  it  follows  a  tortuous  path 
through  the  branches  of  the  tree.  Let  n(x)  represent  the  number  of  pores  of  radius  rp  at 
location  x  in  a  tree  of  trunk  radius  rt.  An  analysis  (Simons,  1982)  of  this  pore  tree  has 
demonstrated  that 


n(x)=r?/r*(x)  (  5) 

and  the  coordinate  x  is  related  to  rp  by 

drjdx^rjl,  (1® 

The  continuous  branching  model  has  been  used  to  successfully  describe  char  oxidation 
(Lewis  and  Simons,  1979;  Simons,  1982  &  1983a),  coal  pyrolysis  (Simons,  1983b  &  1984) 
and  the  catalytic  cracking  of  benzene  by  porous  iron  oxides  (Simons  et  al.,  1986).  It  was  also 
used  to  successfully  describe  sulfur  sorption  (S02  and  H2S)  by  porous  calcine  (CaO)  in  the 
limit  of  zero  utilization  (Simons  and  Rawlins,  1980;  Simons  et  al.,  1984)  and  was  later 
extended  to  include  CaS04  and  CaS  deposits  (Simons  and  Garman,  1986;  Simons  et  al.,  1987; 
Simons,  1988;  Simons  et  al.,  1988).  The  subsequent  determination  of  the  controlling  physical 
parameters  led  to  a  new  concept  for  the  optimization  of  the  sulfur  sorption  process  (Simons, 
1991;  Simons  et  al.,  1992)  through  spray  drying  of  water  soluble  organic  calcium  solutions  to 
control  the  sorbent  pore  structure. 


6 


III.  INTERCONNECTIVITY 


The  first  step  in  determining  the  size  distribution  of  the  interconnected  pores  and  the 

distribution  of  the  permeability  is  to  determine  the  distribution  function  Gt(rrr^drp  which 
represents  the  number  of  pores  of  radius  rp  (within  size  range  drp)  per  unit  cross  section  of  an 
arbitrary  plane  and  also  contained  within  a  tree  whose  trunk  radius  is  rt.  Consider  an  infinite 
homogeneous  isotropic  porous  .media  and  isolate  a  spherical  volume  of  that  media  denoted  by 
the  radius  a.  Such  a  volume  is  illustrated  in  Fig.  2.  The  total  number  of  pores  of  radius  rp 
(within  size  range  drp)  intersecting  plane  AA  of  area  7ta2  has  previously  been  defined  by 

g(r^ita2drp  .  The  pores  in  plane  AA  in  this  size  range  may  also  be  determined  by 

integrating  Gt(rt,rp)na2drp  over  all  trees  whose  trunk  intersects  the  exterior  surface  of  the 
porous  sample.  Hence  it  follows  that 


S(rp)^a2drp  -  f  [Gt(rt , r^a2dr^ g(r) 4na2drt 


(17) 


where  only  those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp . 


A  solution  to  Eq.  (17)  for  Gt(rt,r^  will  not  necessarily  be  unique.  Physical  arguments 

will  help  determine  Gt(rt,r p  and  help  ensure  that  it  is  the  particular  solution  we  seek.  Since 
Nt  represents  the  number  of  pores  of  size  rp  in  the  tree  and  the  probability  of  a  pore 

intersecting  a  plane  is  proportional  to  its  length,  it  follows  that  Gt(rt,r p  should  be 


Figure  2.  Spherical  Volume  of  a  Porous  Media 
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proportional  to  the  product  of  Nt  and  l/it  ,  i.e.,  proportional  to  rt2/rp3.  Eq.  (17)  is  identically 
satisfied  by  a  function  which  differs  from  rt2/rp3  by  ln(rp). 


GMprp)=- 


47ta2rjln(rmM/rp 


(18) 


Note  that  ln^p/r,^)  introduces  an  integrable  singularity  at  rp=rmax  such  that  Gt(rt,rp)drp  is 
finite  at  rp=rmax  .  Hence,  there  is  one  and  only  one  largest  pore  for  each  reference  sphere. 

Further  confirmation  of  the  form  for  Gt(rf,rp)  may  be  made  through  the  definitions  of 
porosity  and  internal  surface  area.  If  the  porosity  and  internal  surface  area  in  plane  AA  due  to 
a  single  tree  of  radius  rt  are  defined  by 


and 


6,  =  /  Gt(rt,rp)2nrpdrp 


‘min 


(19) 


(20) 


respectively,  then  the  total  porosity  and  internal  surface  area  must  be  represented  by  their 
respective  integrals  of  0t  and  s,  over  all  trees.  Hence  0 


and  Sp 


0  =  /  Qtg(rt)4%a2drt 


(21) 


sp  =  f  stg(r)4nazdrt 

raia 


(22) 


follow  respectively. 

It  is  stated  without  further  proof  that  Eqs.  (19)  to  (22)  are  compatible  in  both  their 
exact  and  their  approximate  forms.  The  distribution  function  Gt(rt,rp)  as  given  by  Eq.  (18) 


is  exact  and  will  be  used  to  determine  the  interconnectivity  of  the  pore  structure. 

The  probability  of  trees  sharing  common  branches,  i.e.,  the  interconnectivity  of  the 

pore  structure  is  described  in  Fig.  3.  We  seek  the  distribution  function  I(r^drp  which 
represents  the  number  of  pores  (within  size  range  drp  about  rp)  per  unit  area  of  plane  AA  that 
are  connected  to  both  sides  of  the  pore  structure  through  pores  at  least  as  large  as  rp.  \  is 
defined  as  the  area  within  plane  AA  that  is  open  to  one  side  of  the  porous  media  through  all 
trees  of  size  rt'  (through  all  pores  of  size  rp'  that  are  at  least  as  large  as  rp).  Subsequently, 

AQGt(rt,rpdrp  represents  the  number  of  pores  of  size  rp  (within  size  range  drp)  per  unit  area 
of  plane  AA  that  are  contained  in  a  tree  of  size  range  drt  about  rt  and  are  also  connected  to 
the  opposite  side  of  the  porous  media  through  all  trees  denoted  by  r,\  It  follows  that  the 
distribution  function  for  interconnected  pores  in  plane  AA  may  be  obtained  by  integrating 

AQGt{rt,r^drp  over  all  trees  (r,)  that  are  large  enough  to  contain  a  pore  of  size  rp  .  Hence, 


Tirana2 drp  =  J  \AQGt(rt,rp)drp\g(r)2Tca2drt 
r? 

From  the  above  definition  of  A^  ,  A^  may  be  expressed  as 

rnmx  rt 

Ae=  f[f  2nrp  7ta2G{(rr,rp drj g(r)2iza2drt 


(23) 


(24) 


Figure  3.  Interconnectivity  of  a  Porous  Media 
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where  the  primes  on  the  variables  of  integration  have  been  omitted.  Evaluating  Eq.  (24) 
yields 


^0 


%a2Q  ln(r 


■Jrp) 


2(3 


from  which  Eq.  (23)  yields  the  common  branch  distribution  function. 


(25) 


61 njr^lrj 
4p 


(26) 


It  has  been  deduced  that  the  total  number  of  common  branches  of  size  rp  in  an 
arbitrary  plane  scales  approximately  with  the  total  number  of  pores  of  that  size  in  that  plane. 
Hence,  there  is  a  probability  of  interconnectivity  at  pore  size  rp  that  is  logarithmic  in  pore 
size.  Defining  this  probability  as  Pt(rp)  via  Eq.  (26), 


6  In (r^/rp 

4p 


(27) 


it  is  apparent  that  approximately  one  percent  of  all  pores  of  all  sizes  are  interconnected 
through  larger  pores.  Hence,  convection  will  be  responsible  for  transport  in  approximately 
one  percent  of  the  small  pores.  While  permeability  is  dominated  by  the  largest  pores,  it  is 
important  to  determine  the  level  of  convection  that  is  occurring  in  smaller  pores  in  order  to 
accurately  describe  the  fine  scale  transport  necessary  to  assess  chemical  reactions.  The  broad 
size  range  associated  with  the  interconnectivity  suggests  that  a  very  wide  range  of  pore  sizes 
separates  the  "mobile"  and  "immobile"  regions  of  the  pore  structure  and  that  a  complicated 
mixture  of  convective  and  diffusive  transport  persists  in  this  pore  size  range.  Determining 
permeability  as  a  function  of  pore  size  is  the  first  step  in  describing  the  transport  in  this  pore 
size  range. 


IV.  CONVECTION  IN  THE  PORE  TREE 

Convection  across  plane  AA  in  Fig.  4  will  possess  contributions  from  three  sources 
illustrated  in  Figs.  4a,  4b  and  4c.  Fig.  4a  illustrates  the  case  where  the  convection  in  plane 
AA  is  due  solely  to  the  pores  that  are  interconnected  in  that  plane.  Fig.  4b  illustrates  the  case 
where  the  convection  in  plane  AA  is  due  to  the  smaller  pores  in  the  pore  tree  that  are 
interconnected  outside  of  plane  AA.  This  connectivity  will  translate  into  a  slower  velocity  in 
the  pore  crossing  plane  AA  but  could  be  significant  because  99%  of  the  pores  in  plane  A  A 
are  not  interconnected  in  that  plane.  Fig.  4c  illustrates  the  case  where  the  convection  in  plane 
AA  is  due  to  the  pores  in  the  largest  trees  that  provide  a  direct  link  across  the  porous  sample. 
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b)  Pores  Interconnected  Out  of  Plane  AA 


c)  Pore  Tree  Spans  the  Porous  Sample 
Fig.  4  Convection  in  the  Pore  Tree 
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The  fundamental  relationship  governing  convection  in  a  porous  media  is  Darcy’s  Law 
which  relates  the  volume  flow  rate  Q  to  the  pore  radius  rp  and  the  pressure  gradient 

6-^±[  &)  (28) 

V  8p  {  dx) 

where  p  is  the  viscosity  of  the  fluid.  Considering  an  isolated  pore  tree  with  n  pores  at 
location  x  within  the  tree  and  drp/dx  as  given  by  Eqs.  (15)  and  (16)  respectively,  the  volume 
flow  rate  for  the  entire  tree  (n  branches)  is  given  by 


2  3 

*r,rP 

8  \ilt 


dP 


dr 


pj 


For  a  constant  volume  flow  rate  within  the  tree,  integration  of  Eq.  (28)  yields 


(29) 


P(rt )  ~p(rp)  = 


4&WK/  1 


(30) 


which  indicates  that  extremely  large  pressures  will  be  required  in  the  tree  trunk  to  maintain 
the  constant  volume  flow  rate  in  the  smaller  pores.  This  suggests  that,  whenever  possible,  the 
flow  will  "short  circuit"  through  the  largest  interconnected  pores,  the  pressure  gradient  will 
be  depicted  by  the  largest  pores,  and  there  will  be  a  small  but  finite  flow  in  the  smaller  pores. 

For  the  case  of  rt  »  rp  and  p(rj  »  p(rp)  ,  Eq.  (30)  predicts  that  the  volume  flow  rate  of 
the  tree  is  limited  by  the  smallest  value  of  rp  to  the  second  power 


*  r?P(rt)  r] 

4 


(31) 


while  that  of  a  single  pore  is  limited  by  the  value  of  rp  to  the  fourth  power. 

6  _  Q™  _  *P(Qrt  (32) 

p  n(x)  4  \ilt 

Comparing  Eq.  (32)  to  Eq.  (28)  suggests  that  the  pressure  gradient  is  depicted  solely  by  the 
largest  pores. 
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(33) 


( 


dP  =  ~2P(rt) 
dx  lt 

Hence,  flow  in  the  pore  tree  is  consistent  with  that  for  a  bundle  of  isolated  pores  and  the 
volume  flow  rate  within  a  single  pore  tree  may  be  used  to  determine  the  dominant  transport 
modes  identified  in  Fig.  4. 

Consider  any  pore  of  radius  r3  in  plane  AA  of  Fig.  4b  to  be  the  trunk  of  a  tree.  Each 
pore  of  size  rp  within  the  tree  possesses  the  probability  Pj(rp)  of  being  interconnected  and 

each  interconnected  pore  in  the  tree  will  carry  volume  flow  rate  Qp(.rpyrs) •  Since  there  are 
N9drp  (Eq.  10:  N9=r93/rp4)  pores  in  size  range  drp  within  the  tree,  the  total  volume  flow  rate 

Q«,(rs)  through  trunk  r9  in  plane  AA  becomes 

rs 

OJr.)  =  /  Qp(rp,rj  P,(rp)  N ,  drp  M 


or,  to  first  order, 


QJ'J 


e«p(f,)r>(W/r,)  HQT 
16  P  ls 


(35) 


Within  this  approximation,  it  is  seen  that  Q„(rs)  is  identical  to  the  volume  flowing 
through  the  pores  that  are  interconnected  within  plane  AA.  i.e., 


QJO  -  Qp(rs)  P/rs) 


(36) 


which  demonstrates  that  all  volume  flow  through  plane  AA  in  pore  size  r9  is  dominated  by  the 
interconnectivity  of  size  r9  in  plane  AA  and  not  by  the  interconnectivity  of  smaller  pores  in 
subsequent  branches  of  the  pore  tree.  Simply  stated:  case  4a  dominates  case  4b. 

Case  4c  has  not  yet  been  evaluated. 
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V.  SUMMARY 


This  report  is  the  third  technical  report  on  the  project  entitled  "Development  of  Pore 
Structure  Models  for  Water  and  Contaminant  Transport  in  Partially  Frozen  Soils,"  submitted 
to  ARO  for  the  period  1  November  1994  to  18  November  1994.  The  pore  tree  model  is  being 
extended  to  describe  the  permeable  pore  structure  which  characterizes  the  subsurface  flow  of 
water  in  soil,  the  dispersion  of  contaminants  and  the  in-situ  remediation  of  contaminated  sites. 
Permeability  requires  a  statistical  determination  the  "branches"  that  are  common  to  several 
trees  to  allow  percolation  through  the  large  scale  (mobile)  structure  in  addition  to  diffusion 
through  the  smaller  scale  (immobile)  structure.  While  permeability  is  dominated  by  the  largest 
pores,  it  is  important  to  determine  the  level  of  convection  that  is  occurring  at  the  intermediate 
scales  in  order  to  accurately  describe  transport.  Approximately  one  percent  of  all  pores  are 
interconnected  through  larger  pores.  Hence,  convection  will  be  responsible  for  transport  in 
some  of  the  small  pores  without  influencing  the  over  all  permeability.  This  suggests  a  very 
wide  range  of  pore  sizes  separating  the  mobile  and  immobile  regions  and  a  complicated 
mixture  of  convective  and  diffusive  transport  in  this  pore  size  range.  Subsequent  analysis  will 
utilize  the  interconnectivity  of  the  pores  to  determine  the  distribution  of  the  permeability  with 
pore  size. 
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